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Magnetic Phase Diagram and Metal-Insulator Transition of NiS2 x Se a 
Masato MATSUURA*, Haruhiro Hiraka, Kazuyoshi Yamada** and Yasuo Endoh 

Department of Physics, Tohoku University, Sendai 980-8578 
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Magnetic phase diagram of NiS2- :c Se :E has been reexamined by systematic studies of electrical 
resistivity, uniform magnetic susceptibility and neutron diffraction using single crystals grown by 
a chemical transport method. The electrical resistivity and the uniform magnetic susceptibility 
exhibit the same feature of temperature dependence over a wide Se concentration. A distinct 
first order metal-insulator (M-I) transition accompanied by a volume change was observed only 
in the antiferromagnetic ordered phase for 0.50 < x < 0.59. In this region, the M-I transition 
makes substantial effects to the thermal evolution of staggered moments. In the paramagnetic 
phase, the M-I transition becomes broad; both the electrical resistivity and the uniform magnetic 
susceptibility exhibit a broad maximum around the temperatures on the M-I transition-line 
extrapolated to the paramagnetic phase. 

KEYWORDS: metal-insulator transition, mott-insulator, antiferromagnetism, neutron diffraction, phase diagram 
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§1. Introduction 

N1S2 has been recognized for a long time as a proto- 
type of Mott insulator described by the Mott- Hubbard 
narrow band picture. Hence, many experimental results 
on magnetism and transport properties in this system 
have been interpreted by this scenario, 1 ) where a metal- 
lic phase appears by a broadening of 3d-bands with cither 
applying pressure beyond ^3.5 GPa or ~25% substitu- 
tion of S sites with Se atoms. 

However, a recent photo-emission spectroscopy experi- 
ment revealed that subsequent change in the band struc- 
ture upon Se-substitution should not be described by the 
band broadening but is governed by the charge-transfer 
within a band gap. 2 -* Moreover, recent observations by 
the angle-resolved photoemission spectroscopy (ARPES) 
in the metallic phase near the phase boundary of metal- 
insulator transition (M-I boundary) showed a sharp peak 
in the spectrum near the Fermi energy. 3,4 ) Enhancement 
of the effective mass near the M-I boundary 5 ' also sug- 
gests a breakdown of the previous band picture. 

An antiferromagnetic (AF) long-range order has been 
confirmed to exist in both metallic and insulating 
phases. 6 - 1 The AF long-range ordered structure in N1S2 
is complicated showing two types of structure (c|mi = 
(001) and q M 2 = (323) on fee lattice). 7 " 10 ) Although 
several studies on magnetic phase diagram of NiS2- K Sec; 
were carried out, 6,11-14 ) there still exist non-trivial dis- 
crepancies in the Neel temperature (Tjnj) and M-I tran- 
sition temperature (Tmi)- Furthermore, no systematic 
neutron study on NiS2- a; Se a ; with single crystal has been 
carried out. 

We succeeded in growing a series of single crystals over 
a wide Se-concentration. Then, we performed a system- 
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atic study to reexamine the magnetic phase diagram and 
the M-I transition itself, and established a new phase 
diagram of NiS2-xSea; system with which an intimate 
relation between magnetic and transport properties was 
made it clear. 

§2. Experimental Detail 

Prescribed amount of Ni, S and Se powders in 5N pu- 
rity were mixed with 3% of excess S and Se. They were 
sintered in an evacuated silica tube at 720°C for 7 days. 
In order to improve the homogeneity of starting powders, 
we repeated the above process of sintering and grinding 
twice. The Se-concentration of polycrystalline samples 
was determined from the lattice constant by utilizing Ve- 
gard's law between NiS 2 (5.688A) and NiSSe (5.813A). 6 ) 

Then, single crystals were grown by a chemical trans- 
port method using CI2 gas. 15 ) Polycrystalline powder of 
0.5 g was sealed under 0.5 atm of CI2 gas in an evac- 
uated silica ampoule (100 in inner diameter and 14-20 
cm in length) and was placed in a furnace. The average 
temperature was kept at about 750°C with the temper- 
ature gradient, 2°C/cm. Single crystals were grown up 
to 4 mm on an edge for a month of growth. They are 
characterized by shiny (100) and (111) facets. 

The electrical resistivity (p) was measured by the stan- 
dard four-probe method between 4.2 K and room tem- 
perature. The uniform magnetic susceptibility (x) was 
measured with a standard SQUID magnetometer under 
a magnetic field of 1 T from 5 K to 300 K. Neutron 
diffraction measurements were performed on the triple- 
axis spectrometer TOPAN installed at the JRR-3M Re- 
actor of Japan Atomic Energy Research Institute. Inci- 
dent and final neutron energy was fixed to be 14.7 meV 
(A = 2.67 A- 1 ) using pyrolytic graphite (PG) monochro- 
mator and analyzer. Horizontal collimation of the neu- 
tron beam was set to be 60'-30'-60'-100', from the fore- 
front of the monochromater to the entrance of neutron 
detector. The sample was mounted so as to access to 
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the (h,h,l) reciprocal lattice plane. In this paper, we de- 
note the indices of reflection in the AF unit cell which 
is twice as large as the chemical unit cell, as used in 
previous works. n ) Temperature dependence of lattice 
constant was measured by using a standard four-circle 
spectrometer with x-rays provided from a rotating an- 
ode (40 kVxlOO mcV, MoKai) and monochromized us- 
ing PG(002) reflection. 

In order to investigate the M-I transition in detail, we 
performed simultaneous measurements on p and stag- 
gered magnetization for x = 0.50 and 0.53 by contacting 
copper wires for resistivity measurement onto the crystal 
mounted in an aluminum can for the neutron diffraction 
measurement. Intensities of AF(002) Bragg reflection 
(type I) and p were monitored every 30 and 50 seconds, 
respectively. We changed the temperatures at a constant 
rate of 0.5 K/min. 

§3. Experimental Results 

3.1 Neutron diffraction 

Bragg reflections corresponding to type I AF struc- 
ture were observed in both the insulator and the metallic 
phase (x > 0.5), while reflections of type II AF structure 
vanish for x > 0.3, which is consistent with the previous 
result. 16 ) Thermal evolutions of peak intensity of (002) 
AF Bragg reflection I(T), as typical examples of type 
I, are shown in Figs. 1 (a) for x = 0.50 and (b) for 
x = 0.69. 

6000 r 




(a) 
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0.35~0.40, which is consistent with the theoretical value 
of the three dimensional Heisenberg model (0.367). On 
the other hand, for samples with 0.44<i<0.59, ft exceeds 
0.5. We note that the unusually large (3 is observed only 
for samples with close to the M-I boundary. 

Assuming the noncollincar type I AF structure, 9 ' 10 ) 
we calculated the magnitude of sublattice magnetic mo- 
ment from the intensities of two AF (002), (220) and 
nuclear (222) reflections measured at 10K. We used the 
diffraction data from polycrystalline powder samples to 
minimize the uncertainty of extinction effect. The mag- 
netic moment thus obtained is shown in Fig. 2. Our re- 
sult almost reproduces previous one by Miyadai et al. 1 ^ 
However, our data indicate a change in Se-concentration 
dependence around the M-I boundary; the staggered mo- 
ment decreases more rapidly with Se-substitution in the 
metallic phase than in the insulating one. 




0.2 0.4 0.6 0.8 1 

NiS^Se, x 



Fig. 2. x-dcpcndcncc of the magnetic moment for NiS2— x Se x at 
10 K. Closed and open circles represent present data and those 
by Miyadai et al., 16 ^ respectively. Dotted line shows the M-I 
boundary at K. 



Fig. 1. Thermal variation of peak intensity of (002) antifcrromag- 
netic Bragg reflection for NiS2-^Sea; (a) :r=0.50 and (b) x=0.69 



We determined both and a critical exponent (3 on 
the basis of eqs.(l) and (2), assuming a Gaussian distri- 
bution of Tn . 

7(T)oJ J(t?) 2 '/^)^ (T<T N ) (1) 
I (T>T N ), 

t(rp \ 1 ( (TNO-TN) 2 ^ 

The least-square fitting gives — 45.0 K, a = 7.4 
K and [3 = 0.71 for x = 0.50, and T N0 = 89.4 K, 
a = 0.9 K and (i = 0.39 for x = 0.69. The distribu- 
tion of Tn c(Tn) can be converted into the distribution 
of Se-concentration o~(x) in the sample by using the in- 
dependence of Tn. The o-(x) is determined to be <~ 0.01 
for both x = 0.50 and x = 0.69. For samples with x<0.44 
and x>0.60, the critical exponent (3 is determined to be 



3.2 Electrical resistivity and uniform magnetic suscep- 
tibility 

The temperature dependences of p and \ are shown in 
Figs. 3 and 4, respectively. We note that both measure- 
ments were performed for the identical crystals. From 
thermal evolutions of p and x, we can categorize four 
compositional regions in terms of M-I transition; (1) 
semiconducting region(0 < x < 0.47), (2) first-order M-I 
transition (0.50 < x < 0.59), (3) broad M-I transition 
(x ~ 0.65), and (4) metallic region(x > 0.69). 
(1) Semiconducting region (0 < x < 0.47) 

p exhibits a typical feature of activation type for usual 
semiconductors. The activation energy E a at room tem- 
perature decreases with increasing Se-substitution from 
80 mcV (NiS 2 ) to 50 meV (x = 0.47). These values 
are consistent with the hopping energies determined by 
Kwizera et al. 17 ^ 

X shows a well-defined cusp corresponding to de- 
termined by neutron diffraction measurements. is 
almost constant at 40 K in the semiconducting region. 
For MS2 , a weak ferromagnetism appears below T c = 30 
K. For samples with 0.36 < x < 0.47, T c drops down 
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Fig. 3. Temperature dependences of the electrical resistivity (p) 
for NiS2-^Sea; single crystals. Solid arrows indicate T N de- 
termined by neutron diffraction measurements. Dashed arrows 
point at the temperature where p drops discontinuously. Open 
arrows represent the temperature at which p reaches the maxi- 
mum. 



to 20 K with the weak bulk ferromagnetization, which is 
10~ 3 times smaller than that of N1S2. 

p exhibits a plateau or a broad maximum at ?max as 
indicated by open arrows in Fig. 3. T max is higher than 
7n an d decreases with increasing x from 120 K (N1S2) to 
80 K (x = 0.47). x does not follow a simple Curie- Weiss 
law above T^. These facts suggest that the short-range 
magnetic correlation develops well above Tjj, strongly- 
coupled with transport properties. 
(2) First order M-I transition (0.50 < x < 0.59) 

At high temperatures, p shows an activation type be- 
havior. E a determined at room temperature decreases 
with Se-substitution from 30 meV (x = 0.50) to 20 meV 
(x = 0.59). 

The whole process of the M-I transition in NiS2- a; Se a ; 
seems to occur in a broad temperatures range with a 
width of 20^30 K. Both p and x exhibit a maximum 
(open arrows in Figs. 3 and 4) at temperatures higher 
than Tn (solid arrow), but no cusp appears at in x- 
Below Tn, p rapidly decreases with decreasing temper- 
ature, and then further drops discontinuously by order 
of 10 _1 ~ 10~ 2 , indicating the first order M-I transition 
occurs in the AF ordered state. Tyn is defined as the 
temperature where p discontinuously drops. The differ- 
ence between Tyn and becomes small with increasing 
x. 

p further decreases below Tmi, and eventually exhibits 
a weak temperature-dependent metallic behavior. Wc 
have investigated the detailed feature of the M-I transi- 
tion by simultaneous measurements on the electrical re- 
sistivity and the staggered magnetization, as described 




T(K) 



Fig. 4. Temperature dependences of the uniform magnetic sus- 
ceptibility (x) for NiS2-a;Sea; single crystals. Solid and dashed 
arrows indicate T N and T MI , as in Fig.3. Open arrows represent 
the temperature at which \ reaches the maximum. 



in §3.3. 

(3) Broad M-I transition (x ~ 0.65) 

Although the phase transition from semiconductor to 
metal occurs on cooling, no discontinuous change was 
observed in p and \. As shown for x = 0.65 in Figs. 
3 and 4, both p and x exhibit a broad maximum in the 
paramagnetic phase; at the temperatures higher than Tn 
by about 60 K. T max increases with increasing x, while 
In saturates in this region. In contrast to the region (2), 
X shows a small cusp at T^. p exhibits metallic behavior 
at low temperatures. It is noted that this broad M-I 
transition crossovers to the following metallic one at the 
Se-concentration with the maximum T]\j. 

(4) Metallic region (x > 0.69) 

p shows a typical metallic behavior below 300 K as 
shown in Fig. 3 for x > 0.69. However, a broad max- 
imum in x still exists well above Xn as indicated by 
open arrows in Fig. 4, which suggests the persistence 
of spin correlation at much higher temperatures above 
In. Similar to the previous work, 18 ) T 2 dependence in 
p was observed at low temperatures. Crossover from T 2 
to T-linear behavior occurs around temperatures where 
X shows a maximum. 

3.3 Detailed simultaneous scan on electrical resistivity 
and neutron diffraction near M-I transition 
As described above, the first order M-I transition was 
observed for 0.50<a;<0.59 in the AF ordered phase. In 
order to study correlation between the M-I transition and 
the magnetism in more detail, we performed simultane- 
ous measurements on the electrical resistivity and stag- 
gered magnetization; the peak height of (002) AF Bragg 
reflection i"o02- Temperature dependences of p and /q02 
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measured simultaneously are shown in Figs. 5 for (a) 
x = 0.50 and (b) x = 0.53 single crystals. 

For both x = 0.50 and 0.53, p sharply drops below 
Tj<j and the first order M-I transition occurs in the AF 
ordering state. For x — 0.59 as shown in Fig. 5(c), both 
magnetic order and metallic phase appear almost simul- 
taneously. As is clearly seen in Figs. 5, there exists an 
unusual correlation between the staggered magnetization 
and electrical resistivity, p gradually decreases associ- 
ated with tailing in 7o02 around Tn, which gives unusu- 
ally large critical exponent (3. In addition, I002 show a 
small but clear kink at Ty\± where p drastically change in 
Figs. 5(b) and 5(c). Moreover, there exists an unusual 
kink in /002 at temperature below which p and \ show 
metallic behavior. Below the kink temperature, Tk m k, 
temperature dependence of /002 becomes weak in metal- 
lic phase. This phenomenon is not due to the extinction 
for the single crystal, since the same phenomenon was 
also observed for x — 0.53 in the polycrystalline powder 
sample. 
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Fig. 5. Temperature dependences of the electrical resistivity and 
the peak intensity of (002) antiferromagnetic Bragg reflection, 
T-aaiiT) measured simultaneously for (a) x = 0.50 and (b) 
x = 0.53, and separately for (c) x = 0.59. Circles and trian- 
gles indicate electrical resistivity and Ioo2{T), respectively. 



multi-step like change in p. For the data from x = 0.50, 
the multi step change near 30 K in Fig. 5(a) was con- 
cluded to be ascribed in the following extrinsic reason. 
The single crystals of x = 0.50 were grown from a mix- 
ing of two powder samples synthesized separately, which 
gives rise to a slightly difference in concentration. After 
this simultaneous measurement, we synthesized the poly- 
crystalline powder sample of x = 0.50 and grew again the 
single crystal which shows a single step change at the M- 
I transition. Therefore, we concluded the single crystal 
of x = 0.50 showing two step anomalies contained two 
regions with slightly different Se-concentration by about 
x = 0.005. 

However, for x=0.53 and 0.59 we found no distinct 
experimental reason to induce inhomogeneity of Se- 
concentration. The same multi step change of p was 
also observed in a previous independent study on single 
crystals grown by a flux method. 19 ^ Therefore, the multi- 
step change for x=0.53 and 0.59 can be intrinsic unlike 
the present data for x = 0.50. We only speculate that 
between T^nk and Tmj, the metallic and the insulating 
phases possibly coexist due to the first-order transition 
at Imj. The multi step change of p may relate with such 
an inhomogeneous state in this temperature region. 

As shown in Fig. 6, the lattice constant discontinu- 
ously changes with a small thermal hysteresis only at 
the M-I transition, but no anomaly at Ik m k. The small 
difference in magnitude of lattice constant between Fig. 6 
and its inset is due to the slight shift in position of sam- 
ple, 0.02°, on cooling process. This first order M-I tran- 
sition can be interpreted to be associated with the vol- 
ume contractions with decreasing temperature. Apply- 
ing the Clausius-Clapeyron relation to the present data, 
entropy change at the M-I transition AS is calculated to 
be 1.6 (±0.3) [J/(mol deg.)], as previous report of M-I 
transition induced by pressure. 2 °) This value is in good 
agreement with that from the specific heat measurement 
by Sudo et al, AS = 1.0 [J/(mol deg.)] for x = 0.51 
polycrystalline powder sample. 21 ) 




30 40 50 60 70 



T(K) 

Fig. 6. Temperature dependence of lattice constant for x = 0.53. 
Open and closed circles represent data on heating and cooling 
processes, respectively. 



It is noted p shows a multi step change near the M-I 
transition as shown in Figs. 5. Here we describe on the 
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§4. Discussion 

In Fig. 7, we present a revised phase diagram of 
NiS2-a;Sea; as the summary of the present study. Here, 
we are mainly concerned about the phases near the M- 
I phase boundary for x > 0.3. Compared with pre- 
vious results of the phase diagram, 6 ' n " 14 ) we add two 
crossover regions as shown in Fig. 7. In the lower Se- 
concentration range, the crossover region depicted by 
vertical lines in the figure is determined by the temper- 
ature variation of transport properties. For example, p 
exhibits a broad maximum at the crossover region. Such 
a crossover was also observed in thermoelectric power 
measurements. 17,19 - 1 Furthermore, our recent neutron in- 
elastic scattering data for MS2 show that the short-range 
magnetic correlation remains above at least up to ~ 
200 K. Therefore, we conjecture that the crossover be- 
havior in the transport properties is inherently correlated 
with the magnetic short-range correlation. We call the 
region below the crossover as 'short-range antifcrromag- 
netic insulating (AFI)' phase. The onset of short-range 
AFI phase decreases with increasing x, while is almost 
constant in the insulating phase. 

Another crossover in the higher Se-concentration range 
approximately corresponds to the M-I transition-line ex- 
trapolated towards the paramagnetic phase. Around 
the crossover temperature p and \ show a broad peak. 
For x > 0.69, the peak in p seems to disappear. For 
0.50 < x < 0.59, the peak appears at temperature 
slightly higher than that of the first order M-I transi- 
tion. Here, we add a new phase called anomalous metal 
between Tn and T max . In Fig. 7, we refer the data 
by Ogawa for x > l.O. 22 -* We note that such a broad 
maximum in \ was a ^ so observed in the paramagnetic 
phase of La 2 -£ C Sr 2: Cu04, high temperature supercon- 
ducting cuprates. 23 ) 
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Fig. 7. Magnetic phase diagram of NiSg-^Se^. Open circles rep- 
resent T N . Bold line shows a line of first order M-I transition (sec 
Fig. 8 for more detail). Lines are guide to the eyes. The region 
depicted by vertical lines indicate a crossover region where trans- 
port properties or x(T) reach a maximum; p(T) and Seebeck co- 
efficient 17 ' 19 ) for x < 0.47, p(T) and X (T) for 0.50 < x < 0.65, 
and x{T) f° r x > 0.69. For x > 1.0, we referred the data by 
Ogawa. 22 ) AFI and AFM denote antiferromagnetic insulator and 
metal, respectively. 



A detailed magnetic phase diagram near the M-I 
boundary is shown in Fig. 8. We confirmed that the 
first order M-I transition occurs only in the AF ordered 
phase for 0.50 < x < 0.59 from the simultaneous scan 
on p and neutron diffraction, and the M-I transition be- 
comes second order in the paramagnetic phase, as is seen 
in p and \- According to the recent band calculation 
by the two-band Hubbard model in infinite dimension, a 
first order M-I transition in long range AF ordered phase 
is expected to occur in consideration of the electron hop- 
ping between nearest neighbor Ni 3d and ligand 3p. 24 ^ 
It suggests the charge-transfer nature being indispens- 
able to understand the mechanism of M-I transition in 
this system. We note that the charge-transfer nature is 
seen in ARPES spectrum as a sharp peak near the Fermi 
energy in the metallic phase around the M-I boundary. 

It is important to compare the M-I transition in 
NiS2- a; Se a ; studied here with the previous experiment 
under pressure for MS2. Both phase diagrams under 
pressure and Se-substitution are qualitatively the same 
with each other. 25 > However, there exists an important 
difference between the two. In fact, a recent electrical 
resistivity measurement revealed that the first order M-I 
transition in MS2 occurs at 150 K at 3.0 GPa, 26 ) which 
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Fig. 8. Detailed magnetic phase diagram near M-I boundary. 
Solid squares point at temperatures where p(T) shows a kink. 
Well-defined first order M-I transition was observed at T M i, rep- 
resented by X. 
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can be converted into x ~ 0.6 for Se-substitution. On 
the other hand for Sc-substitution, the maximum tem- 
perature for both the first order M-I transition and the 
long range AF ordering are limited to 100 K. Therefore, 
the observation of the M-I transition at 150 K under 
pressure may conflict with our conjecture that the first 
order M-I transition occurs only in the antifcrromag- 
netic phase. We speculate however, local randomness 
by Se-substitution may reduce the transition tempera- 
tures. Therefore, it is highly required to determine the 
7n under high pressure beyond ^3.5 GPa for N1S2. 

The revised phase diagram in NiS2- a; Se 2: is quite sim- 
ilar to those of other strongly correlated electron sys- 
tems, such as V2O3, 27 ) and organic conductors (BEDT- 
TTF^X, 28-30 ) where antiferromagnetism and M-I tran- 
sition are closely connected. In the latter material, 
a superconducting phase appears, while V2O3 and 
NiS2_ :c Se a ; become AF metals. 

In these systems, it is often observed that the on- 
set temperature of long range ordered state is signifi- 
cantly suppressed compared to the thermal robustness 
of short-range magnetic correlation. This is the case for 
NiS2-^Se x . Since Ni 2+ ions form an fee lattice in this 
structure, the nearest neighbor spins can be magnetically 
frustrated. In fact, preliminary inelastic neutron scatter- 
ing measurement on this system found an unusual mag- 
netic excitations possibly due to this frustration effect. 

§5. Summary 

NiS2- a; Se :E system has been reinvestigated systemat- 
ically by neutron diffraction, electrical resistivity, uni- 
form magnetic susceptibility and X-ray diffraction using 
both single crystals and powder samples. We revealed 
a clear correlation between the transport and magnetic 
properties in this system. Both the magnitude of mag- 
netic moment and its temperature dependence are af- 
fected by electric transport properties, particularly by 
the M-I transition. A well-defined first order M-I tran- 
sition was found to occur only in the antiferromagnctic 
state. On the other hand, the transition in the paramag- 
netic phase becomes broad and is seen in the resistivity 
up to the Se-concentration where Tn exhibits a maxi- 
mum. The novel feature in the newly proposed magnetic 
phase diagram is the existence of two cross-over regions, 
one from short ranged AF phase to paramagnetic insula- 
tor (PI) and the other from anomalous metal to PI. These 
crossover features suggest that the short-range magnetic 
correlation remains well above Tn in both phases, which 
cannot be interpreted by a single Mott Hubbard model, 
but by a new model taking account of the strong electron 
correlations. 
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